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ABSTRACT 



Context. T-type dwarfs present a broad and shallow absorption feature centred around 6950 A in the blue wing of the K doublet 
^ ' at 0.77 fim which resembles in depth and shape the satellite absorption predicted by detailed coUisional broadening profiles. In our 

^SJ I previous work, the position of the predicted line satellite was however somewhat too blue compared to the observed feature. 

On ■ Aims. In this paper, we investigate whether new calculations of the energy surfaces of the potentials in the K-H2 system, including 

spin-orbit coupling, result in a closer coincidence of the satellite with the observed position. We also investigate the extent to which 
CaH absorption bands contribute to the feature and at what T^ff these respective opacity sources predominate. 

Methods. We present model atmospheres and synthetic spectra, including gravitational settling for an improved description of depth- 
dependent abundances of refractory elements, and based on new K-H2 line profiles using improved interaction potentials. 
Results. By comparison with a high signal-to-noise optical spectrum of the Tl dwarf elndiBa, we find that these new models do 
reproduce the observed feature, while CaH does not contribute for the atmospheric parameters considered. We also find that CaH is 
settled out so deep into the atmosphere that even turbulent vertical mixing would appear insufficient to bring significant amounts of 
CaH to the observable photosphere in dwarfs of later type than ~ L5. 
. ■ Conclusions. We conclude that previous identification of the feature at this location in the spectra of T dwarfs as well as the latest L 

' dwarfs with CaH was erroneous, as expected on physical grounds: calcium has already condensed onto grains in early L dwarfs and 

?H ', thus should have settled out of the photosphere in cooler brown dwarfs. This finding revokes one of the observational verifications for 

the cloud-clearing theory assumption: a gradual clearing of the cloud cover in early T dwarfs. 

Key words. Stars: low mass, brown dwarfs - Stars: atmospheres - Line: profiles 



1. Introduction 

Late L- and T-type brown dwarfs have atmospheres composed 
primarily of molecular hydrogen and helium. Most refractory 
metals condense out to grains in previous evolutionary phases 
(early L type) and should settle below their fully radiative up- 
per photospheres. Alkali elements bind less readily to molecules 
or grains and their resonance transitions remain the last sources 
of opacity at optical wavelengths, along with Rayleigh scatter- 
ing by H2 and He. Alkali line broadening therefore defines the 
local "continuum" out to several thousands of Angstroms from 
the line cores of the n eutral K and Na D doublets a t 0.59 and 
0.77 jum, respectively. iBurrows & VolobuvevI (I2003L and refer- 
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ences therein) have shown that these far wings deviate greatly 
from the simple Lorentzian shape, and speculated on the exis- 
tence of satellite features in the blue wings due to a local ex- 
tremum in the energy shift. 

A number of T dwarfs sh ow an absorption feature at 6950 A 
which lBurgasser et alJ (120031) identified as a CaH band. This was 
a somewhat curious suggestion, however, because calcium con- 
denses onto grains in L dwarfs already and is thus expected to 
have settled out of the photosphere of cooler brown dwarfs. The 
identification of this feature with CaH (as well as the observation 
of FeH bands) led to believe that it was a proof for vertical mix- 
ing and clearing of the cloud cover in T dwarfs ( iBurgasser et all 
,2003). 

As alternative, we propose that the feature should rather 
be identified with quasi-molecular absorption by K-H2. 
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Fig. 1. FORS2 red optical spectrum of the Tl dwarf eIndiBa (solid), compared to our synthetic spectra for a 1.3 Gyr (dashed) and 
2 Gyr (dotted) model. To highlight the various sources of opacity, spectra obtained when the K and Na D doublets are omitted (dot- 
dashed), and when dust grains involving Ca and other refractory species are prevented from forming and raining out (upper dotted 
line) are also shown. All molecular bands but CaH and H2O have been omitted from the latter spectrum for clarity. The EFOSC 
spectrum of eIndiBa,b is shown from 8600- 10000 A (light [magenta] solid line) against a composite model of both T dwarfs. 



In. Allard et alj (l2003bl) calculated absorption profiles of the 
K doublet at 0.77 /im for H2 and He perturbers, based upon 
the K-H2 and K-H e interaction poten tials calculate d by 
iRossi & Pascal^ (I1985L hereafter£P850 and lPascald ( Il983[ here- 
afterlP83h using a semi-classical approach. The resulting models 
and synthetic spectra revealed a K-H2 satellite feature with sim- 
ilar depth and width to that observed at 6950 A in the T dwarfs, 
but at a shorter wavelength of 6830 A. In this paper, we present 
new models based o n fine profiles for K-Hg calculate d wi th up- 
dated potentials (see lN. Allard et al.Ll2007al hereafter [S06() . 



2. Model Atmospheres 

We have used th e stell ar atmosphere code PHOENIX 
dHauschildt & Baronl 1 19991) version 15. With respect to 
previous versions used bv I Allard et al] (1200 Ih to determine the 
limiting effects of dust grain formation of the atmospheres of 
low-mass stars and brown dwarfs, the chemical equilibrium 
(CE) code in th is version has been m odified to include grain 
settling effects dF. Allard et al.L l2007cl) . For each layer of the 
model atmosphere, atarting at the deepest cloud condensation 
layer (as set by the CE calculation), dust grain number densities 
are calculated in equilibrium to the gas phase. Condensation, 
gravitational settling and turbulence mixing timescale are 
com pared to predic t the fra ction of grains which must settle (see 
also iF. Allard etall l2003i . This grain fraction is then removed 
from the composition and a new equilibrium is obtained. This 
process is repeated until the grain density no longer changes. 
The advantage of this method compared to an a priori reduction 
of refractory elemental abundances is an automatic depletion 
of the elemental abundances as a function of temperature and 
pressure, in the sequence of formation of the diverse grain types. 
For example, Zr will be depleted first from the gas phase, as 
Zr02 is the grain that forms at the highest temperature, followed 



at lower temperature s by Fe , Mg (MgSiOs), Al (AI2O3), and so 
on. The Allard et al.l (1200 lb "Cond" models — which used an 
equilibrium dust distribution but simulated gravitational settling 
by commenting out dust opacities - tended to retain some trace 
of the molecular compounds involved in the grains. Here the 
depletion is complete and, for example, no TiO absorption band 
opacity remains in layers where grains involving Ti have been 
found to condense. This formalism corresponds therefore to a 
more accurate description of grain settling. 

In the Settl models described above, gravitational settling ef- 
ficiently clears the uppermost region of the photosphere from 
dust down to a few pressure scale heights from the top of the 
convection zone. A full description of these models exceeds the 
scope of this paper and will be published separately. 

3. New line profiles 

An important improvement with respect to the 2001 Cond mod- 
els is the inclusion of the pressure br oadening profiles f or the 
neutral alkali resonance fines of Li i bv All ard et £11(12005 ), Rb I 
and Cs i bv .All ard & Spiegelma n (2006), K i and Nai D as in 
In. Allard et aLr (l2003b) and bv In. Allard et all (l2007al) . These 
profiles include fine structure and use damping constants for the 
Lorent z line cores from semi-classical theory dN. Allard et al.L 
I2007bh . 

The most important contributions to the satellite feature un- 
der consideration come from H2 and He. The profiles of K-H2 
account for two different orientations of the H2 molecule: C21., 
where the H2 molecular axis is perpendicular to the K radiat- 
ing atom, and Ccov, where all atoms are collinear The C21. IS06l 
satellite is predicted to lie at 6980 A, while the weaker Ccov con- 
tribution peaks at 6920 A: these should be compared to 685 1 and 
6695 A, respectively, for the IRP85I profiles. The K-He satelHte 
is predicted to be at 6930 A. Note, however, that when account- 
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ing for intermediate orientations of the H2 molecule, additional 
satellite contributions may wash out and slightly change the pre- 
dicted position. S antra & Kirby (2005) have published ab ini- 
tio potentials for K-H2 at 4 orientations of the H2 molecule, as 
well as for K-He, predicting K-H2 satellites around 6920 A and 
6850 A for the C21. and Ccav orientations, respectively, with K-He 
at 7000 A. These values correspond to the extrema of the differ- 
ence potentials, however: our corresponding values are 6906 A 
and 6 830 A for C21. and Ccov, respectively. Finally, IZhu et all 
(l2006h have also computed profiles of K and Na perturbed by 
He, finding a K-He satellite at 7080 A. Combined, these calcula- 
tions predict a broad absorption feature with properties very sim- 
ilar to that observed at 6850 to 7 100 A and centred at 6950 A) in 
the spectra of T dwarfs such as s Indi Ba as discussed below. 

4. Optical spectral analysis of the T1 dwarf £ Indi Ba 

The eIndiBa,b binary system is a unique test for brown dwarf 
models. They are the closest known brown dwarfs to the Sun, 
allowing spectra to be obtained with high signal-to-noise in the 
faint optical range. The well-determined HIPPARCOS distance 
(3.626 pc) allows accurate comparison to be made with models 
at an ab solute level. Based on their JHK photometry and //-band 
spectra, McCaughrean et al. (2004) obtained effective tempera- 
tures of 1276 K and 854 K for Ba and Bb, respectively, using the 
iBaraffe et al.l ( l2003l) 1 .3 Gyr Cond isochrone for the most likely 
age of the eIndi AB system. This corresponds to surface gravi- 
ties of logg = 5.15 and 4.9, and radii of 0.091 and 0.096 Rq, re- 
spectively. We have compared a Settl model calculated for these 
parameters to a medium-resolution resolved optical spectrum of 
the Tl primary (Ba), as shown in Fig.[T] 

As part of a comprehensive 0.6-5/im photometric and spec- 
troscopic study of s Indi Ba,b, the ESO VLT was used on June 16 
2004 to obtain optical (0.6-0. 86;um) spectroscopy using FORS2 
in long-slit mode, with a slit width of 0.5 arcsec and the HR 
colUmator. In this mode, the 600RI grism yielded a spectral res- 
olution of R~1000 at 6780 A, while the seeing of 0.35 arcsec 
FWHM ensured that the two T dwarfs, separated by ~0.8 arcsec 
at that epoch, were well-resolved spatially. A total integration 
time of 80 min was obtained by co-adding 6 x 800 s individ- 
ual exposures. The data reduction was standard, using skylines 
for wavelength calibration, dome lights for spectral flat-fielding, 
and the DC white dwarf spectropho tome trie standard, LTT 9491, 
for flux calibration. Unfortunately, no direct telluric calibration 
was possible, but an excellent cancellation of all known tel- 
luric features was achieved by using the N SO Kitt Peak atmo- 
spheric transmission spectrum ( Hinkle et al.ll2003h convolved to 
the FORS2 resolution and with iterative adjustment of the effec- 
tive airmass to match the conditions over Paranal at the time of 
our observations. This allowed efficient removal of the 6884 A 
O2 B band close to the satellite (compare Fig.|2]l. Further details 
of the data reduction will be given by King et al. (2007, in prep). 
Preliminary absolute flux calibration was obtained by using the 
combined /-band p hotometry from the DENIS second data re- 
lease Q, as listed in iMcCau ghrean et all ( 120041 ). to calibrate a 
composite spectrum of sIndiBa.b taken with the EFOSC spec- 
trograph at the ESO 3.6 m telescope, the red part of which is 
also shown in Fig. [1] The FORS2 spectrum was then rescaled 
to match the composite spectrum at A < 8000 A, where the flux 
contribution of the Bb component is less than 1 %. 
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Fig. 2. Effective temperature sequence of logg = 5.0 Settl mod- 
els with (solid) and without (dashed) CaH band opacity. The syn- 
thetic spectra have been degraded to 5 A resolution for compar- 
ison with L dwarf spectra ( thick Tred] lines) observed with the 
Keck LRIS spectrograph ( Reid et all 1200 iL 12000; Marti n et all 
1999). Spectral subtypes o f the L dwarfs are given in the opti- 
cal classification system of iKirkpatrick et all (1 19991 in prepO); 
no te fliat DENIS J 025 5-4700 has been typed L6 in the system 
of lMartmetaLl (|1999|). The sharp CaH band head near 7075 A 
vanishes in both observed and model spectra between about L5 
and L8, and a distinctive absorption feature emerges that extends 
well to the red of the original CaH band head. The models have 
Li abundance reduced to 1 % (solid) and 0.1% (dashed), indicat- 
ing strong Li depletion in all these field L dwarfs. 
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Our synthetic spectrum reproduces the near-IR photometry 
for Ba to within 0. 15 mag, supporting our treatment of dust opac- 
ity, which shapes the infrared SED in transition objects. The red 
optical region is extremely well reproduced in most details, in- 
cluding the shape of the K doublet at 0.77 /im and the depth of 
the secondary alkali lines of Rb i and Cs i, absorbing through the 
red wing of the K i doublet. Assuming a meteoritic, i. e. proto- 
solar, lithium abundance, the model also predicts a strong Li i 
6707 A doublet, but this is not observed. This implies a rela- 
tively massive object which depleted its Li at a young age when 
its atmosphere was fully convective and its core temperature ex- 
ceeded the Li burning threshold: this process consumes most Li 
in brown dwarfs more massive than 0.055 - 0.06 Mp in the first 
Gyr of their Hfetime (IZapatero Osorio et~aill2005l) . 

According to the iBaraffe et al.l (l2003h evolutionary tracks, 
a 0.055 Mo brown dwarf at an age of 2.0 Gyr should have 
Teif = 1350 K, log g = 5.30 and R = 0.086 Rq. This model is over- 
plotted as dotted line in Fig. [T] and still agrees with the data 
within observational and modelling uncertainties. For a qualita- 
tive estimate of the degree of lithium depletion we have reduced 
the elemental Li abundance in this latter model by 1.5 dex. It 
still shows a clearly visible resonance line not detectable in the 
observation, confirming very low Li abundance. 

The models also show a weak 6894.5 A Sri line against the 
background of the alkali lines. It is not detectable in the observed 
spectrum, indicating incomplete treatment of Sr condensation. 
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Fig. 3. Partial pressures of principal Ca-bearing species versus 
optical depth in the Settl model for T^^ff - 1280 K (thick lines) 
and 1600 K (thin lines). Gas temperature is shown (dotted lines) 
with the convection zones indicated (crosses), right scale. 



4.1. The K-Ho satellite 

The key feature in this high signal-to-noise spectrum for present 
purposes is the shallow absorption feature centred at 6950 A, the 
exact location of the satellite of the K-H2 interaction with the 
0.77 jum doublet as newly predicted here using the |S06| profiles 
(Figs. [Hand 12]). However, the new satellite apparently extends 
not far enough to the red to reproduce the more elongated shape 
of the observed feature. K-He produces a satellite which might 
reproduce the observed ex tent of the feature to 7 100 A accord- 
ing to the lSantra & Kirbvl (|2005) and lZhuetall (|2006|) calcula- 
tions. Our old potentials predict it further to the blue, from 6930- 
6980 A, thus superimposing on the location of the K-H2 satellite 
and producing an absorption trough that is too narrow. 

4.2. The CaHA^n-X^I. band 

The dot-dashed line in Fig.[T]shows the eff'ect of eliminating the 
large wings of the K and Na D doublet, to reveal other opacity 
sources which participate in the formation of the spectrum. The 
Settl model does not predict any features besides the Sri line 
(compare also Fig.|2]l in the spectral region of the satellite. 

However, if we artificially prevent calcium from being 
locked into grains and settled out, the A^n - X^l, band of CaH 
is seen, partly coincident wi th the observed f eature (upper dot- 
ted line in Fig. [TJ. This led iBurgasser et alj (l2003h to propose 
that CaH is responsible for the observed feature in early to mid- 
T dwarfs, advocating a possible resurgence of CaH due to up- 
welling of CaH or cloud-clearing effects. 

We find that most calcium is locked into diopside 
(CaMgSiaOe) and akermanite (CaaMgSiiO?) below 1900 K, 
leading to a depletion of calcium to 10""* or less as these conden- 
sates settle out (Fig. [3]l. CaH only prevails at solar abundances 
deep in the optically-thick atmosphere, independent of whether 
we account for clouds or not. 

Fig. |2] illustrates the difference between models with and 
without CaH bands for effective temperatures ranging across the 
late-L sequence. At 2000 K, the CaH band does indeed shape 
the pseudo-continuum at 6800-7050 A, but by 1600 K, the K-H2 
satellite is dominant and CaH barely contributes to the absorp- 
tion. While the intensity of the satellite actually increases with 
gas temperature (see Figs. 10-12 of S06), in our model spectra 
the satellite grows in strength with decreasing Teff , due to clear- 
ing of the atmosphere and increased pressure in the line-forming 



region (see IF. Allard et al.L l2007d for details). It creates an ab- 
sorption trough extending from 6850 to 7 100 A, visible in the 
observed spectra of the late L dwarfs, where it is extending be- 
yond the CaH band head and even further to the red than in the 
models (cf. § 14.11 ). Once T dwarf temperatures are reached, our 
models retain essentially no CaH opacity, challenging the cloud- 
clearing picture unless condensable gas is upwelled efficiently 
over more than two pressure scale-heights by turbulence. 

5. Summary 

We have computed new model atmospheres based on detailed 
absorption profil es f or the neutral NaD, Li, K, Rb, and Cs al- 
kali lines (|RP85| an d |P83l potentials) and, in particular, new in- 
teraction potentials iSOo) for K-H2. These new models predict 
a K-H2 satellite absorption feature at 6950 A closely matching 
the position and shape of an observed feature in the spectrum 
of the Tl dwarf eIndiBa. We therefore conclude that the K-H2 
satellite is the most natural explanation for the feature, rather 
than CaH as has been previously proposed. Indeed, the high (so- 
lar) abundance of CaH which would be required to generate ab- 
sorption near this spectral location should only be found deep in 
the optically-thick atmosphere at the reff=1280K of eIndiBa, 
preventing its bands from appearing in spectra even assuming a 
clearing of the surface cloud coverage. 
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